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The strong dependence of the momentum distribution of the photoelectrons on experimental con- 
ditions raises the question as to whether angle-resolved photoemission spectroscopy (ARPES) is 
able to provide an accurate reflection of the Fermi surface in Bi-based cuprate superconductors. 
In this paper we experimentally prove that the main contribution to the intensity variation comes 
from matrix elements efi^ects and develop an approach to overcome this problem. We introduce a 
concept of 'self-normalization' which makes the spectra essentially independent of both the matrix 
elements and particular experimental parameters. On the basis of this concept we suggest a simple 
and precise method of Fermi surface determination in quasi-2D systems. 
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I. INTRODUCTION 

Since the beginning of the field of the high Tc supercon- 
ductors (HTSC), ARPES has taken a_apecial role in the 
experimental study of these systems.EjD Among the nu- 
merous remarkable ARPES experiments on the HTSG, a 
special place belongs to the investigation of the Fermi sur- 
face (FS) of these systems. The difficulties encountered 
upon applying 'traditional' techniques for determining 
the FS to the HTSC (such as de Haas - van Alphen and 
positron annihilation) focussed the attention on the pos- 
sibility offered by ARPES to obtain a direct image of the 
basal-plane projection of. the FS. The vast majority of 
the earlier ARPES workoiZi agreed in the view that the 
FS in the most HTSC is hole-like, and centred at the X,Y 
points of the 2D Brillouin zone. Recently, however, a con- 
troversy regarding the FS topology of these systems has 
flared up. Some groups have suggested the presence of 
an electron-like, F-centred FS in the Bi2Sr2CaCu208+5 
(BSCCO) [Refs. 1-0] and Pb-BSCCO [Ref. |ll| systems, 
which would represent a complete revision of our thinking 
regarding the fundamentals of the electronic structure of 
the HTSC. On the OjthetJiand, other groups have also re- 
visited this questionO^Ej and confirmed the 'old' picture 
of a hole-like FS in the Bi-cuprates. At the same time 



in La2-a;Srj;Cu04 a cross-over from hole-like to electron- 
like FS has been suggested from ARPES data_p|B_going 
to the overdoped side of the phase diagramme.Ejil3 Con- 
sequently, there still exists no consensus as to the correct 
picture for the normal state Fermi surface topology in the 
HTSC, making this question an important one to clarify. 

The current debate as regards the Fermi surface topol- 
ogy i&^|baS|e:L.to a large extent upon ARPES intensity 
maps.ElEHijilj This means that the issue of matrix ele- 
ments, which could be strongly photon enei:fl';^^nd k- 
dependent in the 2D cuprate-based materialsE£l~Lj has to 
be treated seriously. Thus we should be able to identify 
situations in which the matrix elements dominate and, 
where possible, develop practical methods of exctracting 
the underlying true information from the raw photoe- 
mission intensity signal. A second issue is that of how 
to accurately determine the Fermi momentum vectors, 
ki? from real photoemission data. The accuracy within 
which this k^ determination can be tested has increased 
dramatically in the past few years as a result of a new 
generation of electron energy analysers which offer res- 
olutions in k-space one order of magnitude superior to 
what was previously available. 

In this paper, we address the question as to how one 
can best locate Fermi momentum vectors in the HTSC 
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with the aid of the angular distribution of photoemis- 
sion intensity. As a case study we take BSCCO, but in 
fact the conclusions arrived at are quite general to the 
high resolution ARPES investigation of quasi-2D sys- 
tems. We experimentally demonstrate the strong dis- 
tortion of the 'pure' photoclcctron angular distribution 
caused by matrix element effects, thus making an appro- 
priate further analysis of the raw data compulsory. To 
give such a statement a firm foundation, a detailed dis- 
cussion of the experimental conditions and their influence 
on the photocurrent together with a critical overview 
of the existing methods of hp determination are pre- 
sented. We then demonstrate that, using what we call 
a 'self-normalization' procedure, one can significantly re- 
duce the dependence of photoemission spectra on the ma- 
trix elements and finally show that this approach can be 
successfully applied to the BSCCO compounds. Conse- 
quently on the basis of the self-normalization method we 
formulate a criterion of determining the Fermi surface of 
the HTSC from ARPES data. 



II. METHODOLOGY 

Because of rapid evolution of the modern ARPES ex- 
periment and rising number of possible techniques for the 
solution of a given problem, we include in this section not 
only traditional experimental subsection (A) but also a 
quite detailed description of our approach to the FS map- 
ping. In subsection B we discuss the quantities which are 
in principle accessible by ARPES in our implementation. 
Then we discuss how to optimize the experimental pa- 
rameters for the study of the HTSC cuprates (subsection 

C) before closing this section with a critical evaluation of 
the available methods of kj? determination (subsection 

D) . 

A. Experimental 

Two types of experimental set-up have been used. The 
majority of the data discussed here were recorded with 
an overall energy resolution of 19 meV (FWHM) using a 
SCIENTA SES200 analyser coupled to a high intensity 
He resonance source (GAMMADATA VUV5000) via a 
toroidal grating monochromator (giving a degree of linear 
polarization of ca. 40%). The SES200 analyser provides 
an angular resolution of down to 0.2°. The single crystals 
were mounted on a purpose-built, high precision cryo- 
manipulator which allows the sample to be rotated with 
a precision of better than 0.2° about three perpendicular 
axes in a wide range of angles. The synchrotron based 
data were recorded as described in Rcf. 



B. ARPES with analysis of the energy and 
momentum distributions on an equal footing 

The new generation of electron energy analysers men- 
tioned above have enabled a jump in angular resolu- 
tion performance as a result of electron optics possess- 
ing an angular dispersion capabilitv-ija the direction par- 
allel to the analyser entrance slit.E3 One can visualise 
the new mode of data collection with the help of Fig. 
|l|. This shows the information landing on the analyser's 
2D-detector within 7 minutes' measuring time in a static 
experiment - i.e. without moving the sample or the anal- 
yser. One direction on the detector represents an angular 
interval and the other direction - an energy interval (in 
this case ±7° and ~0.6eV respectively). 

The dataset shown in Fig. |l| is a 'snapshot' taken 
along the r-(7r, tt) high symmetry direction of the Bril- 
louin zone of Pb-doped BSCCO at 120 K, with the pho- 
toemission intensity plotted as a function of both bind- 
ing energy and momentum. The left panel shows cuts 
of the intensity distribution /(k, lu) parallel to the en- 
ergy axis, i.e. energy distribution curves (EDC), which 
are uniquely defined by the fixed value of momentum. 
One can also cut the same /(k, w) distribution parallel to 
the momentum axis. These cuts are shown in the right 
panel and are termed momentum distribution curves or 
MDCs.EI An MDC should reflect the vector nature of the 
momentum, and thus is uniquely defined by the chosen 
frequency (binding energy) and an arbitrary path in a 
two-dimensional k|| -space. It is clear from Fig. that 
the modern ARPES machinery enables the simultane- 
ous measurement of the energy and angular distribution 
of the photoelectrons leaving the sample and therefore 
permits a treatment of the energy and momentu p]. . d| e- 
pendence of the photocurrent on an equal footing.aO 

The use of a grey scale to represent the intensity en- 
ables an efficient presentation of the three dimensional 
data set shown in Fig. |l|. The result, shown in the up- 
per panel of Fig. |^ is called an energy distribution map 
(EDM). An EDM can be thought of as an array of EDCs 
(or MDCs) taken along particular path in the BZ within 
particular range of binding energies. If we now fix bind- 
ing energy and grey-scale-code the intensity in a series of 
MDCs covering an area in (k^,, kj,)-space together, then 
we arrive at a momentum distribution map or MDM, 
which represents a constant energy surface. The lower 
panel of Fig. |2| shows such an MDM for E—0 eV binding 
energy covering a part of the first Brillouin zone of Pb- 
BSCCO. In our case, the sample rotation involved in the 
recording of an MDM is such that the individual MDCs 
representing the angular breadth of the 2D detector are 
radial in nature, with the origin at the F point. 

The informative capacity provided by such a map is 
easy to estimate, although not only the Ep-MDM is im- 
portant. The dataset still possesses the binding energy 
axis, and so a series of MDMs then represents the evo- 
lution of the momentum distribution of the electronic 
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states when going-from the Fermi energy towards higher 
binding energies Ej 

The inset to Fig. ^ summarizes the completeness of the 
information available in our ARPES experiment. Here we 
illustrate the three dimensional (k^:, ky, ti;)-space which 
can be probed with high E and k resolution. The fourth 
dimension here is symbolised by the grey scale and repre- 
sents the photoemission intensity. Recording the ARPES 
intensity while moving alcpg any vertical direction, i.e. 
parallel to the energy axisEj will give an EDC, whereas 
an MDC is the intensity distribution along the arbitrary 
path which belongs to any of horizontal planes in this 
space. Horizontal planes themselves are MDMs and the 
vertical surface (plane) defined by a given path (line) in 
an MDM is the EDM. Thus^ the portion of (k, a;)-space 
shown in the inset to Fig. is confined by three EDMs 
(one is not visible) and two MDMs (one is not visible) and 
consists typically of approximately 100000 data points. 

The foregoing discussion has illustrated the potential 
offered by our experimental setup. Nevertheless, despite 
significant advances there remain factors such as the res- 
olution or lifetime of the sample surface which make it 
necessary to optimize the other experimental conditions 
for the treatment of the physical problem at hand. More- 
over, given the widespread use of intensity maps in the 
literature regarding the FS topology of the HTSC, it is 
evident that the chosen experimental conditions (such 
as the experimental geometry, the excitation energy, the 
photon polarization, the temperature etc.) could deci- 
sively change the final picture obtained. Therefore, in 
the next sub-section we briefly deal with the different 
experimental parameters which could strongly influence 
the photoemission intensity distribution in (k, u;)-space. 

C. Factors influencing the /(k,(j;) distribution of 
photoelectrons 

In order to discuss different parameters affecting the 
measured photoemission intensity, we first write an ex- 
pression for the photocurrent as a function of (k, lu)- 
space. For quasi-2D systems and under ±he assumptions 
that the 'sudden approximation' applie£3 and that only 
a single initial state is involved, the photocurrent can be 
written in the following form 

/(k,c^) = Gk{M{k)[A{k,Lu)f{Lo)] ® i?^,k -t- B{lo)}, (1) 

where Gk is a mainly geometrical prefactor which will 
be described below, M represents the square of the ma- 
trix element linking the initial and final states, A is the 
single particle spectral function, / is the Fermi function 
and Ru,,k is the energy and momentum resolution func- 
tion. B is the background, which contains extrinsic ef- 
fects such as inelastic scattering of the photoelectrons. As 
an approximation we assume a negligible k-dependence 
of the extrinsic background and likewise a negligible uj- 
dependence of both the matrix elements and prefactor 



Gk within the energy interval of interest (~ 0.3 eV). Eq. 
1 makes it clear in a formal manner, that the measured 
signal is not simply the spectral function, and thus that 
a number of parameters must be known before A can be 
extracted. 

The pre-factor Gk describes the combined effects of 
extrinsic parameters such as occur upon the rotation of 
the sample with respect to the analyser (i.e. changing 
effective photon density in the area of the sample 'seen' 
by the analyser) or the inequal efficiency of the different 
channels of the parallel detection system. The raw-data 
MDM shown in Fig. ^ (lower panel) illustrates the ef- 
fect of Gk, as it can be seen that at the interjoins of the 
two separate arcs of individual (radially arranged) MDCs 
there is an intensity misfit, mainly due to inequal detec- 
tor channel efficiencies. In order to minimize the effects 
connected with Gk, calibration scans can be carried out 
by measuring an isotropic photoemitter such as an amor- 
phous gold film. Another, simpler way of overcoming this 
problem, is the self-normlaization which will be described 
later. 

The generic step-like background B observed in the 
ARPES of the HTSC is still a puzzle. The authors of 
Ref. ^ demonstrated that the contribution from sec- 
ondary electrons, which could be estimated wpon the ba- 
sis of electron energy-loss spectroscopic datao, is not suf- 
ficient to explain the background which is in-correspon- 
dence with earlier assumptions in this regard.EHI However, 
our assumption that the background is approximately 
k-independent and of practically negligible intensity at 
Ep in comparison with the main signal is supported by 
the similar B{uj) line shape for k-points from unoccupied 
part of the BZ or for those k for which spectral func- 
tion peaks at higher binding energies (e.g. close to the 
F-point). Therefore, where needed, the background can 
be safely subtracted. We have found that as a good rep- 
resentative for the B{uj), an EDC from the vicinity of 
{tt/2, tt/2) or (tt, tt) points could be taken. 

The most important component of Eq. |^ other than 
the spectral function itself is the Af(k) term describing 
the matrix elements, which depends upon both the pho- 
ton energy and the photoelectron momentum via the op- 
erator which couples the final and initial state wave func- 
tions. The choice of the energy of exciting photons is 
far from being unimportant even in the case of quasi-2D 
electronic systems. Firstly, upon changing the photon 
energy one alters the momentum resolution. Secondly, 
the 2D-Cu02 plane materials exhibit extremely strong 
variations in the ARPES intensity of their lowest hKriK 
ionisation states as a function of the photon energyc3'EJ 
Thus, in extreme cases, by an unfortunate choice of the 
photon energy the contribution to the angular distribu- 
tion of the photoelectrons from a particular initial state 
can be significantly suppressed. 

The same goes for the 'angular' part of the M(k) 
when using highly polarized radiation. If the experimen- 
tal geometry can be controlled so as to give a clearly 
defined symmetry condition - such as can be the case 
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along high symmetry Hnes in k-space - the strong polar- 
ization dependence of the photoemission signal can be 
used as a_pijobe of the symmetry of the initial states 
involved. EiJEil If, however, the strongly polarised radia- 
tion is used to measure ARPES spectra away from the 
high symmetry directions in the Brillouin zone - then 
the observed photoemission intensity represents only a 
part of the whole picture. In this context we note that 
it has been clearly demonstrated that ARPES intensity 
maps recorded from BSCCO using the same photon en- 
ergy differ very significapt|ly when recorded with differing 
polarization geometries .E3 

Two ways around this problem spring to mind. Firstly, 
one could use unpolarized light. The laboratory He 
source and monochromator used here generates VUV 
radiation with ca. 40% linear polarization. Thus, the 
majority of the MDM intensity is coming from excita- 
tion with unpolarised light, meaning that although the 
polarized component will favour emission from partic- 
ular states, the global effect is quite small and we are 
consequently able to 'see' all the states involved. The 
second possibility is to use the variable polarization of- 
fered by modern insertion devices at synchrotron radia- 
tion sources to record intensity maps in pairs with com- 
plementary polarization geometry. 

Moving further through the factors separating a real 
ARPES experiment from v4(k, uj) we come to the energy 
and momentum resolutions, represented in Eq. |l| by the 
function Ru:.k- When using an angle- multiplexing anal- 
yser, momentum resolution can be projected onto two 
mutually perpendicular directions - parallel to the en- 
trance slit of the analyser and perpendicular to the slit: 
-Ra;,k = -Rw-Rfc||3Ht-^fcisHt ■ The resolution parallel to the 
slit Rk^i^iit is defined by the electron optical character- 
istics of the spectrometer whereas Rk^M is futher con- 
trolled by the aperture and entrance slit size. In most 
cases we used 19 meV x0.2° x 0.5° FWHM resolution. 

Apart from the angular resolution of the analyser, the 
flatness of the sample surface as well as the excitation 
energy and the absolute values of the momenta define 
the momentum component of the (k, w)-resolution. For 
example, the k-resolution for Fermi level emission in ex- 
periments using 'high' photon energies (e.g. 55 eV) is up 
to a factor 3 worse than while using typical 'low' photon 
energies (e.g. 21.2 eV as here) for k-vectors in the second 
Brillouin zone of BSCCO for the same angular resolution. 

The practical effect of the energy and momentum res- 
olutions on the measured data depends strongly on the 
dispersion of the feature in (k, a;)-space. The influence 
of each component increases as the direction of the most 
rapid change of intensity /(k, uj) approaches the corre- 
sponding axis onto which the total resolution is pro- 
jected. The direction of the most rapid intensity varia- 
tion roughly coincides with the normal to the bare band 
(which is a surface in (k, u;)-space). Fig. ^ illustrates 
three exemplary EDCs recorded from pure BSCCO at 
40K. In each case at least one of the resolution contri- 
butions is zero for our experimental geometry. Fig. ||(a) 



shows the kp-EDC taken from the FX cut, for which 
Rk±^iit is negligible - the width of the feature is defined 
jointly by Rui and Rk^^^^-^- In fact, for the strongly dis- 
persing states along the nodal line (~ 2 eVA), the mo- 
mentum resolution along the slit is the dominating factor. 
For the data of Fig. ^(a) the energy resolution was 19 
meV FWHM and Rk^^^,^, was 0.015 A'^ (i.e. 0.2°). The 
latter causes an energetic broadening of some 30 meV, 
which easily outweights the contribution from R^ itself. 
In choosing an EDC from the MX cut as shown in Fig. 
^(b), we switch off the influence of Rk^^,^^^- As in the FX 
case, the momentum resolution (this time Rk^siu) still 
plays the leading role as regards the instrumental broad- 
ening of the observed peak. Fig. ||(c) illustrates an EDC 
from the M-point, where it is well-known that the sharp 
peak observed below, r3jc varies little in binding energy 
as a function of k,Ba meaning that the instrumental 
contribution to the width of this structure is determined 
effectively by the energy resolution. The three examples 
shown in Fig. || show that the conditions we apply result 
in a good balance between all three resolution compo- 
nents, when considered in the light of the typical Fermi 
velocities in the BSCCO-based materials. 

Having discussed the influence of all the parameters 
entering into Eq. |l|, we conclude that the most unpre- 
dictable and therefore difficult factor to deal with which 
separates the photocurrent from the spectral function are 
the matrix elements, which, as was mentioned above, are 
strong in quasi-2D Cu02-plane materials. 

We now turn our attention to a discussion of the best 
manner in which iip vectors can be derived from the 
ARPES data of the HTSC. This point is much more than 
merely a detail of the ARPES data evaluation, and lies 
at the heart of the current Fermi surface controversy. In 
particular, given the conclusion that it is only the ma- 
trix elements that severely hamper a clear view of the 
spectral function, our discussion of how to determine k^ 
is thus centered on the question as to which level each 
method is immune, if at all, from matrix element effects. 

D. What is the best way to determine when k=kF? 

1. Dispersion method (EDC maximum) 

The relatively coarse k-mesh available in the majority 
of earlier ARPES investigations meant that FS crossings 
could only be located by the analysis of a series of EDCs 
(i.e. an EDM) containing a dispersive feature. The most 
simple and intuitively direct method is then to follow the 
energy position of the EDC maximum and to extrapolate 
obtained dispersion relation to E = Ep. This procedure, 
however, suffers from a number of drawbacks. Firstly, the 
influence of the Fermi cut-off distorts the picture within 
ca. 2kT of Ep. Secondly, if the self-energy is frequency 
dependent, following the EDC maxima will lead to the 
wrong result. This point can be simply visualised with 
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the help of the three-dimensional plot, shown in the cen- 
tral panel of Fig. |l|. As the intensity varies along the 
bare band, the trace formed by joining the maxima of 
cuts through this object taken parallel to the binding 
energy axis (EDCs) can never agree with that obtained 
by joining the maxima of the cuts taken parallel to the 
momentum axis (MDCs). If the self-energy only weakly 
depends on k, it is evident that the MDC dispersion is 
much closer to the 'true' dispersion (i.e. the bare one plus 
the real part of the self-energy) than that from the EDC 
maxima. As a consequence, even though it does possess 
the advantages that it is insensitive to the normalization 
procedure, to the effects of finite energy resolution and is 
quite robust with respect to matrix elements effects, the 
EDC method gives only approximate values of kp. 



2. A r method 

It has been proposed fropa, ARPES measurements of 
TiTc2 and from simulationsEEl that Ep-MDCs shift as a 
function of temperature in such a way that the differ- 
ence of such MDCs turns out to be zero only for k=ki?. 
The practical application of this method to the HTSC 
is blocked by two points. Firstly, everything has to be 
measured twice (for Ti and T2) with very high k-space 
location precision, with all other parameters being kept 
equal. This is often impossible due to the finite lifetime 
of the cleaved surfaces of the HTSC. Secondly, and more 
fundamentally, the AT method cannot function if the 
width of the Ef-MDCs concerned is temperature depen- 
dent. As this is very clearly the case in the IITSC,o the 
AT method is invalid in the context of the HTSC. With- 
out wishing to persue this point further here, we refer the 
reader to the Appendix for the analytical evidence which 
forms the basis of these statements. 



3. Symmetrization 

The symmetrization method is based upon an analysis 
of the lineshape of the result obtained by mirroring the 
photoemission EDCs aromidK^Ep and summing up the 
two spectra for each k-pointj23 a procedure that can be 
described as IsCk^u) = /(k, a;) + /(k, — cj). Within this 
method k^? is defined as the point at which the dip at Ej? 
in the symmetrized EDCs (k<kF) turns to maximum. 
It is evident, however, that upon approaching hp the two 
peaks originated from the lu and -uj spectra will approach 
each other and become indistinguishable, giving a maxi- 
mum in the symmetrized EDC over a range of k. One can 
estimate this k-range quantitatively by solving the equa- 
tion {cPls{'k,uj)/duj'^)^=o = 0. For the model spectral 
fimction considered in the Appendix, the solution gives 
k-ki? ^ 0.02 A^^ at 300K. Consequently, the k point at 
which the dip in the symmetrized EDCs transforms into 



a peak is shifted away from the true hp introducing con- 
siderably larger error than, for instance, the MDC max- 
imum method gives (see below). Symmetrization does 
possess the advantage that it eliminates the Fermi cut- 
off from the k^-EDC and answers the question whether 
a FS crossing occured or not in a given EDM even with 
the presence of strong matrix element effects. 



4- Maximal gradient of the integrated intensity 

It is well known that even for interacting Fermi sys- 
tems, hp is characterised by a jump in the momentum 
distribution, n(k). For finite temperatures one could still, 
in principle, detect rapid variations in n(k) and estimate. 
hp from the max |Vk?^(k)|. It has been proposecE3H 
that the integrated intensity {lint ) of an EDC could give 
a measure of n(k) at one particular k-point and thus the 
analysis of a series of EDCs could present an opportunity 
to estimate hp. This method has been applied to differ- 
ent Mfstftms—jncluding BSCCO, and is still intensively 

usediem 

There exist, however, the following arguments against 
the |V/i„t(k)| method. Firstly, /mt(k) is not equal to 
n(k), for the same reasons that the raw photoemission 
intensity is not equal to the spectral function. Secondly, 
a single EDC does not represent the photoemission in- 
tensity fpc, a single k-point, but rather for a range of 
k-points.c3 For lower photon energies (e.g. 21.2 eV), the 
finite k-interval associated with the finite energy width 
of an EDC can even be comparable with the momentum 
resolution. Thirdly, the band structure of the system has 
to be amenable, in the sense that a single state needs to 
be isolated and well away in frequency from other fea- 
tures. Even if this is the case, the intensity integration 
should be carried out over all frequencies - in practice 
however, energy windows varying between 100-600 meV 
in width are taken for the integration. Taking a narrower 
window reduces the similarity with n(k), whereas a wider 
window results in enlargement of the k-interva£3 and in- 
creases the contribution from deeper lying valence band 
states. A further difficulty arises from the data analysis 
in that numerical differentiation introduces additional er- 
rors. Such transformation of the raw data also produces a 
set of additional " false" features on the mapeJ which have 
to be identified as such and neglected at a later stage. In 
any case, the quantitative precision of the hp determi- 
nation is in direct relation to the width of |V/mt(k)|. 
This width is much broader than, for example, a typ- 
ical Ep-MDC for the 2 eVA-like dispersive features in 
the BSCCO compounds. The factors mentioned above 
make the |V/mi(k)| method intrinsically inaccurate. In 
the Appendix, we use simulations to show further that 
this method can result in substantial systematic errors 
determining kj?, .-which is in agreement with the results 
of other authors.Ea 
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5. Maximum intensity method (MDC maximum) 

The maximum intensity method, as first introduced for 
BSCCO in Ref. | is based on measuring the Ep-MDM. In 
this case, the photocurrent is recorded only in a narrow 
energy window centred on the chemical potential, thus 
enabling the coverage of large areas in momentum space 
within a relatively short time. 

The physical basis of this method oi \ip determina- 
tion is straightforward. One starts with the reasonable 
assumption that at finite temperatures the spectral func- 
tion (independent of the model used to describe it) has 
a peak at lo=0 only for k=ki?. It then immediately fol- 
lows from Eq. |^ that every Ep-MDG corresponding to a 
path in an £^i;'-MDM which crosses the Fermi surface will 
show a maximum. This property of Ep-MDCs has been 
recognized and successfull|Y apralicd for \<ip determination 
by a number of groups JIjoSEJ For the MDC peak to 
lie exactly at k^?, the influence of the CJ^-factor, the ma- 
trix elements (M(k)) and resolution (i?w,k) should not be 
strong enough to shift the peak position of the spectral 
function. In this context, the essentially symmetric res- 
olution functions can certainly be regarded as harmless. 

6. Influence of matrix element effects on determination 

It is much more tricky to evaluate how strong the de- 
pendence on the matrix elements is for a particular set of 
experimental conditions. Numerical calculations of the 
photoemission intensity including the matrix elements 
have been carried out, and predict that the matrix el- 
ements should have a dramatic effect on the angular dis- 
tribution of photoelectrons that would be detected in an 
ARPES experimpit, given an identical underlying Fermi 
surface topologyJlj 

In Fig. |(a) we show Ep-MDCs and in Fig. |(b) 
the integrated ARPES intensity lint both recorded for 
a path in k-space along the F-M-Z direction in BSCCO 
compounds. In each case the four panels show data mea- 
sured using different photon energiesji^ whereby the 21 
eV panel was recorded using radiation from a He reso- 
nance source. 

Comparing the shape of the MDCs and lint traces at 
different excitation energies it is clear that a k^ determi- 
nation method which involves the photoemission inten- 
sity such as the maximal gradient of integrated intensity 
or the MDG-maximum methods may function poorly in 
such a caseMj The vertical dashed lines marking the ap- 
parent location of Fermi vectors show different results be- 
tween the two methods, and, more importantly, different 
results for the same method, depending on the experi- 
mental conditions. This holds even for the maximum- 
MDC method. 

This strong dependence of the apparent location of 
ki? on the experimental conditions used to measure the 



ARPES data is, in fact, at the root of the current con- 
troversAjrjeegacdkig the Fermi surface topology of the 
HTSC.iy'EirMH At this stage one could even be led to 
doubt the value of ARPES as a method of determining 
the Fermi surface in the HTSC. Obviously, there is an 
urgent need to find a method which is able to accurately 
refiect the Fermi surface even in the presence of strong 
matrix element effects. In the following, final results sec- 
tion, we describe an approach which allows not only an 
estimation of the distorsions caused by matrix elements 
but also enables a robust and precise determination of 
\<ip vectors. 



III. FERMI SURFACE MAPPING 

From the discussion of the existing methods of k^? 
determination it follows that, in the absence of strong 
matrix element effects, the most precise, simple and 
physically transparent approach is the MDC maximum 
method. Therefore one natural way forward is to improve 
this method by minimizing its sensitivity to the matrix 
elements. The form of Eq. [| suggests the possibility of 
being able to 'divide out' both the Gk and M(k) prefac- 
tors, providing the denominator is proportional to their 
product. Within the energy range under consideration, 
one can take both M(k) and Gk to be frequency inde- 
pendent. Thus, a perfect candidate for such as division 
would be a signal from the same k-point, i.e. from the 
same EDC. At the same time, the reference signal should 
be a slowly varying function of momentum in the vicinity 
of the expected MDM (MDC) maxima. In principle, we 
should restrict ourselves to a rather small energy interval, 
so as give a narrow ki?-window for each EDC. 

Although several possibilities exist for such as matrix- 
element elimination, in the following we concentrate on 
the division by the integrated intensity, as this has al- 
ready been successfully applied to the determination of 
the FS topology in BSCCO systems.N 

As discussed above in the context of the |V/i„f(k)| 
method, the integrated intensity versus k proves to be 
quite a slowly varying function in BSCCO in the vicin- 
ity of the expected Fermi surface crossings in comparison 
with the k-dependence of the Fermi energy intensity (i.e. 
the Ep-MDM). Within the framework of Eq. |^, the in- 
tensity after division, Inorm is given by: 

_ /(k,0) _ 

~ X/(k,c.)dc^ ^ (^j 

^ Gk{A/(k)[A(k,c^)/(a;)] j?^,k|o +-B(0)} 
Gk{/, A/(k)[A(k,L^)/(t^)] (g) R^,^du + I S(u;)du;} ' 

where e is the energy window of integration, normally 
chosen between 600 meV and -100 meV. It is easy to see 
that prefactor Gk cancels out immediately without any 
additional assumptions, thus automatically solving the 
problem of the detector efficiency calibration. 
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As this method involving the division by the integcated 
intensity has been criticised as being 'unphysical'Jlj we 
now consider the behaviour of the function described in 
Eq. 2 above in different parts of the BriUouin zone in 
detail. We consider three regions in the zone: 'definitely 
occupied' (E^ > oJmax, where E'k is the quasiparticle 
dispersion and LOmax is the higher binding energy of the 
integration window), 'definitely unoccupied' {E\^ < 0) 
and 'close-to-the-Fermi-surface' (i?k ~ 0) regions, and 
discuss the contribution to signal at Ep from the spec- 
tral function and extrinsic background. The first region 
contains k-values for which the spectral function peaks 
more than 0.5 eV from the chemical potential. As re- 
gards the 'definitely unoccupied' region, the Fermi cutoff 
means that the signal at Ep from the spectral function 
disappears faster going into the unoccupied part of the 
Brillouin zone, and thus we consider the 'definitely un- 
occupied' part as being located immediately after the FS 
crossing. In both of these regions, the Inorm signal is de- 
termined only by the background lineshape. The weak k- 
dependence of the background means that Inorm should 
be approximately constant in the 'definitely unoccupied' 
and 'definitely occupied' regions, although it is somewhat 
noisy due to the small values of B(0). In the immediate 
vicinity of the expected Fermi surface, we can neglect 
both terms connected with the background, as -B(O) is 
typically an order of magnitude lower than I(ki?,_Bi?) (as 
can be estimated from, for example, the data of Fig. |^) 
and the corresponding ratio of the integrated intensities 
is ^ 0.15. In this case, the Inorm function can be rewrit- 
ten as: 




_B(o) — ^ const 

\ B{uj)duj 

\A(k,uj)f{uj)](g,R^_^\o 
I' lA{k,uj)f{uj)]®R^_kdui 



if Ek ^ 0. 



These relations show that, in the vicinity of k^?, the 
function is independent of both the matrix elements and 
the Gk-factor. It should be noted, however, that in the 
case of severe suppression of the intensity related to the 
spectral function, the background contribution becomes 
substantial and therefore Inorm cannot be considered as 
matrix element independent - even in a first approxima- 
tion - but still can be used for identifying such a situation. 

Although not immediately obvious, the physical mean- 
ing of the Inorm function is quite transparent. Consider 
its behaviour along a single cut through the 2D BZ (e.g. 
that shown in Fig. The numerator is simply an Ep- 
MDC and the denominator is, in this case, the lint along 
FX. We expect the integrated intensity to show a slow 
drop in the region of kj?. At the same k-region the true, 
underlying MDC has a sharp maximum exactly at k^?. 
Thus, the division of the narrow Lorenzian-like MDC 
function by the slowly falling lint function does not even 
result in a significant shift of the maximum, and makes 
itself felt only in the asymmetric shape of the renormal- 
ized MDC. In other words, by normalizing the MDC in 
such a manner, we do not change the position of its max- 



imum, thus preserving its most important property as 
an indicator of k^?. We note that the asymmetry of the 
normalised MDC peak can act as an additional guide in 
determining whether the FS has been crossed from the 
occupied part to the unoccupied or in the opposite direc- 
tion. 

Having now discussed the normalised MDC approach 
on paper, as it were, we now put it to the test. As was 
evident in Fig. |^, the determination of whether there is 
a hp vector along the F-M-Z direction in pure Bi2212 
represents a quite a challenge. In Fig. ^ we show the 
result of dividing each of the panels of Fig. |^(a) with the 
corresponding panel of Fig. ^(b) - i.e. Fig. || contains 
the renormalised £'i;'-MDCs obtained by division of the 
raw MDCs by the lint curves. 

The major advance here over Fig. ^ is that now the 
normalized MDCs are all very similar, which, of course, 
is logical as there can be only one spectral function for 
this system. The nearly perfect coincidence between data 
recorded with a wide range of excitation energies is re- 
markable, and leaves no doubt that it was the matrix ele- 
ments which had led to the differences in the raw MDCs 
seen in Fig. ^. This example shows that in this way 
we can get an estimate of the influence of the most un- 
predictable part of the k-dependent matrix elements - 
i.e. that part which is not due to the symmetry selection 
rules. 

Comparison of the raw [Fig. §(a)] and normalized [Fig. 
H MDCs indicates that the strongest photon energy de- 
pendent variations in the raw intensity, /(k), are in the 
vicinity of the M-point. For hv = 21.2eV, Af(k) is asym- 
metric with respect to M-point, whereby emission from 
the states in the first BZ is favoured. The same asym- 
^^try holds for hiyj-^ 32eV, which is in agreement with 
other experiments.E3 The most noticeable difference for 
32 eV photons is the strong suppression of the spectral 
weight around the M-point itself - i.e. the raw MDC 
has a minimum at M. There is less pronounced M-point 
suppression for hi^ — AOeV, and for hv = bOeV the sup- 
pression is no longer observed and the asymmetry has 
changed to favour emission from the 2nd BZ. This last 
point is a forerunner of the severe matrix elements ef- 
fects at hv = 55eT^ whece-jntensity from the second zone 
is completely dominant.M Thus, as regards the spectral 
weight suppression, the k-interval corresponding to that 
part of the MDC which is 'eaten away' by matrix ele- 
ment effects moves from the second BZ to the first when 
the excitation energy is varied from 21 to 55 eV, thus re- 
sulting in highest intensity suppresion at the M point for 
hv ca. 32 e\jT-|as had been predicted from photoemission 
calculations .lla 

The remaining differences between the self-normalized 
MDCs presented in Fig. || could be due to a number of 
factors. For example, the M-point spectral weight sup- 
pression can be drastic enough (see, for example. Fig. 1(a) 
in Ref. ^ for hv ~ 32ey) to mean that one is no longer 
able to distinguish between the real signal and the back- 
ground. Three of the curves presented in Fig. || are from 
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pristine Bi2212, whose i^&MDM is highly complex in 
the vicinity of the M point .Ej Thus some of the observed 
fine structure couldpbe due to the 'FS crossing' of the two 
diffraction replicas. Ed 

Here we wish to stress that using the self-normalization 
method we are able to get an impression of the effects of 
the matrix elements. To determine the matrix element 
itself in a rigorous way, we would need to know the spec- 
tral function. In this sense, the self-normalized spectra 
shown in Fig. || are not equal to A(k, uj) across the whole 
range of kw-space but are very close to the spectral func- 
tion for k near hp. It is this property of self-normalized 
MDC's which makes them so suitable for FS mapping. 

As a final comment to Fig. || we note that the uni- 
versality of the lineshape of the normalized MDC's along 
FMZ in the BSCCO-based HTSC - in which no sharp 
peaks characteristic of FS crossings occur either side of 
the M-point - strongly suggests the absence of main band 
FS crossings in this direction in the BZ and thereby lends 
weight to the hole-like FS side in the Fermi surface con- 
troversy. 

The foregoing discussion was limited to the self- 
normalization procedure using the integrated intensity 
as the normalising quantity. In Fig. |^ we show both 
raw and normalized MDMs recorded at room tempera- 
ture within 3 hours after the cleavage of a Pb-BSCCO 
(OD 72K) single crystal. 

Firstly, dealing with Fig. ||(b), the self-normalized 
MDM confirms the behaviour expected for the function 
given in Eq. ^. Both the 'definitely occupied' and 'defi- 
nitely unoccupied' regions have approximately the same 
intensity (mid to dark grey [red] tone). On approaching 
the Fermi surface from the occupied side, the signal is re- 
duced (giving dark areas) as here the contribution from 
the signal at the Ep is still small but the integrated in- 
tensity is already quite large. There are two futher possi- 
bilities: if a band crosses the Fermi level (e.g. along F-X) 
a sharp increase in intensity is observed giving a bright 
feature on the map.c3 Alternatively, in the case of the 
F-M-Z cut (see Fig. |), the MDC has a plateau near the 
M-point, reflecting the behaviour of the flat band which 
does not cross the FS but approaches near enough to it 
to contribute to the signal at Ep. This leads to a fairly 
uniform intensity (mid grey [red] tone) around the M- 
point. Upon careful analysis of the locations of the MDC 
maxima in both raw and self-normalised maps, we find 
that there is no detectable shift between the two datasets, 
which confirms again that the lint function varies much 
more slowly than the MDC does in the vicinity of the FS. 
We now turn to Fig. ^(c) , showing the self-normalization 
result using the intensity at high binding energies (in this 
case 250 nieV below E^?) as the denominator. The self- 
nomalised MDM in Fig. o(c) displays all the characteris- 
tics of the map shown in^(b), indicating firstly that the 
high binding energy signal is also sensitive enough to the 
matrix element effects to enable their elimination for k 
near to hp. Secondly, this points to the soundness of the 
assumption that the matrix element effects are insensi- 



tive to energy on the range of 0.5 eV. Thirdly, the good 
agreement between the self-normalization based upon the 
integrated intensity or high binding energy denominators 
proves the physically sound basis of the former procedure, 
in contrast to what is claimed in Rcf. Finally, we re- 
turn to the BSCCO FS controversy, and point out that 
Fig. ^ shows without a doubt that the hole-like Fermi 
surface topology, which can be the only conclusion upon 
looking at our FS maps, is not a product of the integrated 
intensity normalization procedure, but is a robust result. 

To summarise this section regarding the use of self- 
normalization to reduce the strong matrix element ef- 
fects in ARPES of the HTSC we can say: (i) self- 
normalization conserves all the advantages of maximum 
intensity method of k^? determination and comes close 
to the ideal of a method which delivers robust, precise 
results even in the presence of strong matrix element ef- 
fects (ii) the denominator used in the self-normalization 
should be a signal which 'feels' the matrix element effects, 
but which varies relatively slowly in k-spacecZl. We have 
shown both the integrated intensity and the intensity at 
higher binding energies to be two candidates which func- 
tion well, (iii) the effectiveness of this method means that 
one can overcome the doubts raised earlier as regards the 
power of ARPES to determine the Fermi surface: the self- 
normalized MDM represents directly the Fermi surface 
map. 

The self-normalization method also implies a formal 
criterion for the FS determination: a k-point in the 2D 
BZ belongs to the Fermi surface when and only when all 
possible normalized MDCs in its vicinity, except may be 
one (to account for inequal intensity distribution along 
the FS itself), have a local maximum at that point. To 
take into account the finite resolution of the experiment 
one would have need to introduce a tolerance angle within 
which the MDCs could be considered as belonging to the 
FS. 

As regards the extension of the self-normalization 
method to quasi-2D systems other than BSCCO, we see 
no limitations in terms of its validity since the method 
is based on very general grounds. Even in the case of 
very low Fermi velocities the self-normalization effec- 
tively compensates already observable shift of the Ep- 
MDC due to the finite energy resolution, i.e. gives exact 
locations of kj^-vectors. 

IV. SUMMARY 

In this paper we have discussed the factors which sepa- 
rate the data of a real ARPES experiment from the spec- 
tral function, which is a highly topical subject in the light 
of the current controversy regarding the ARPES-derived 
Fermi surface topology in the HTSC. Based on high qual- 
ity ARPES data of the BSCCO and Pb-BSCCO sys- 
tems recorded under a variety of experimental conditions 
(photon energies, degrees of polarization), we have sug- 
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gested a simple method which enables an estimation of 
the strength of the matrix elements effects in the ARPES 
of the HTSC, which at the same time allows the precise 
determination of the Fermi surface even when the matrix 
elements are strongly k-dependent. In this approach, a 
self-normalization effectively immunises the momentum 
distribution curves - whose maxima deliver precisely the 
Fermi wavevectors - against the matrix elements and ex- 
trinsic factors separating the photoemission signal from 
the spectral function for k near to kj?. Consequently, 
the self-normalized momentum distribution map of the 
photoemission intensity at Ep gives the most faithful 
reproduction of the underlying Fermi surface topology 
achievable from real ARPES data, and thus provides easy 
access to the quantitative analysis of the Fermi surface 
in these materials. 



the self-energy we use the following approximation: 
T."{uj,T) = v^(au;)2 + {(3TY with a = 1 and /3 = 2 
and T in energy units), which, as can be seen from Fig. 
0, gives a resonable fit to the experimental data. 

Fig. shows typical experimental Pb-BSCCO FX 
EDM's (for 30 K and 300 K) as contour plots (left pan- 
els), together with the results of the simulation (right 
panels). The quasiparticle dispersion includes the ef- 
fect of the real part of self-energy Yj'{lo,T), but in the 
region of interest near to the Fermi level we consider 
£fc — vpk, where vp is simply the renormalized Fermi 
velocity at cj = 0. We took vp— 2 eVA from the exper- 
imental data. Having described the basis of our simula- 
tions, in the following we analyse how accurate the 'kp 
determined by the different methods is. 
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VI. APPENDIX 



A. Maximum MDC method 

To evaluate the precision of this method quantitatively 
we simulate the FX Ep-MDG according to (|) and (|). 
The results of this simulation are presented in Fig. ^j. 

The first observation is that the error in determining 
ki? is very small using this method. For example, for 
room temperature the error is less than O.OOlA"^ for an 
energy resolution of 19 meV as was used in the experi- 
ment. Even for a resolution of 50 meV, the error is maxi- 
mally O.OOTA"^ Secondly, the shift of the 'observed' k_F 
from the true value is only weakly temperature depen- 
dent, which therefore cannot be considered as an obstacle 
to the use of the 'maximum MDC' method. 



In this section we communicate in a little more detail 
certain points as regards the question of how best to de- 
termine k^. To enable a more quantitative analysis of the 
strengths and weaknesses of the various k^? methods, we 
have generated a simulated dataset, based upon a fit to a 
real FX EDM dataset from BSCCO. Starting from equa- 
tion (1) given earlier and assuming that the problems 
of matrix elements, detector calibration and background 
have been adequately dealt with, we adopt a model in 
which the photocurrent can be calculated as 

I{k,uj)<^[A'{k,uj,Rk)f{uj)]®R^. (4) 

To speed up the calculations, we combine the spectral 
function with the momentum resolution Rk in 

The absence of a strong asymmetry in the FX MDC's 
results in a straightforward influence of the momen- 
tum resolution, which ip-.in contrast to the influence 
of the energy resolution,E3 which is taken into account 
via convolution with the resolution function Ri^(lo) = 
{Ru,y/n)~^exp{—uj'^/R^). For the imaginary part of 



B. AT method 

The original AT proposaS is based upon there being 
a temperature dependence of the position of the MDC 
maxima. Thus, the applicability or otherwise of the AT 
method to the HTSC can also be judged from Fig. ^. The 
first point is that, as discussed above, the T-dependcnt 
shifts of the MDC maxima are vftf small in BSCCO, 
in contrast to the case in TiTe2H3 Furthermore, Fig. 
^(a) also shows clearly that there is no common cross- 
ing point on the right flank of the Ep-MDC^s, which is 
a result of the temperature dependence of the width of 
the MDC's, thus making the AT method inapplicable. 
Finally, even if one assumes a temperature-independent 
width of the Ep-MDC, the accuracy of the AT method, 
6k, is related to the uncertainty in the determination 
of the relative intensities of each of the MDC pairs, SI: 
6k ^ SI/(dI/dk)k=kF- In our case, to reach an accuracy 
of 3 X 10"^ A~^, the SI would have to be less than 3% 
of /, which is beyond most present experimental capabil- 
ities. 



9 



C. Gradient Iint{k) 

Fig. ||(a) shows the resuhs of the shxiulation as regards 
gradient n(k), in which 



^ D.L. Feng, W.J. Zheng, K.M. Shen, D.H. Lu, F. Ronning, 
J.-I. Shi moyama, K. KishiOj G. Gu, D. Van der Marel, Z.- 



dljntjk) 
dk 



dA' 
Hk 



® R^{Lj)duj 



(6) 



for ujmin — —0.1 eV, LUmax — 0.6 eV is plotted for four 
different temperatures. None of the maxima dlint{k) / dk 
he on the Ak=0 hue, indicating a systematic error in the 
determination of kj^. Fig. |^(b) shows the temperature 
dependence of this shift away from the true kj?, A/c, plot- 
ted for different Rk = (10, 30, 60, 100) meV which are 
equivalent to (3, 10, 20, 33)xl0-3 A"! or (0.09, 0.27, 
0.54, 0.90)° of angular resolution. Since our currently 
best intrumental angular resolution is 0.2°, we discuss 
the curve for an angular contribution of 30 meV. Here 
the error at low temperatures is between 0.002 and 0.003 
A"-'^, which is as good as the maximum MDC method 
at these temperatures. For higher temperatures (e.g. for 
T>T* in the HTSC for which the Fermi surface is not 
gapped), the error from gradient Iint{k) has risen to 0.008 
A""'^, some 8 times higher than the corresponding value 
for the maximum-MDC method. 
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''^ The FS maps shown in Fig. ^ show the FS arc around the 
nodal direction with unprecedented clarity. This fact is due 
to (i) an intensity enhancement coming from the symme- 
try selection rules, as in this case the E-vector of the ca. 
40% linearly polarised component of the He I radiation was 
aligned perpendicular to the FX direction and (ii) the high 
overall angular resolution for measurement in this direction 
leading to the large signal to background ratio seen in Fig. 
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the normalization signal was a 'scaled down version' of the 
original _Ef-MDM, as could easily be the case taking the 
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FIG. 1. Middle panel: photocurrent versus binding en- 
ergy and momentum along the r-(7r,7r) direction at 120 K in 
Pb-BSCCO (UD85K). Left panel: energy distribution curves 
(EDC) as parallel intensity profiles corresponding to fixed mo- 
mentum values. Right panel: momentum distribution curves 
(MDC) as intensity profiles corresponding to fixed energy val- 



FIG. 2. Upper panel: energy distribution map (EDM) 
from the r-(7r,7r) direction in the Brillouin zone of Pb-doped 
BSCCO recorded at room temperature. Lower panel: mo- 
mentum distribution map (MDM) of Pb-doped BSCCO, 
recorded at room temperature (raw data). White horizontal 
dashed line represents a kf-EDC, vertical ones correspond to 
an i?_F-MDC. In both cases the grey scale represents the pho- 
toemission intensity as indicated. The inset shows the three 
dimensional (k^,, ky, aj)-space which is probed in ARPES of 
quasi-2D systems. 



FIG. 3. EDCs measured at 40 K from the single cleave of 
pure BSCCO (UD 89K) at different points in BZ illustrating 
the contributions from different components of the resolution. 



FIG. 4. (a) i^F-MDCs (circles) and (b) Ii„t (squares) (from 
500 meV to -100 meV binding energies) from ARPES data 
recorded along the FMZ direction in pristine BSCCO for dif- 
ferent excitation energies: top right, 32 eV; bottom left, 40 
eV; bottom right, 50 eV. The top left panel in each case shows 
analogous data for Pb-doped BSCCO measured using 21.2 eV 
photons from a He lamp. The apparent kf locations from (a) 
the MDC-maximum method and (b) the maximum gradient 
of the integrated intensity are marked with vertical dashed 
lines. 



FIG. 5. Normalized Ep-MBCs from ARPES data recorded 
along the FMZ direction in pristine BSCCO for different exci- 
tation energies: top right, 32 eV; bottom left, 40 eV; bottom 
right, 50 eV. The top left panel shows analogous data for 
Pb-doped BSCCO measured using 21.2 eV photons from a 
He source. 
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FIG. 6. MDM's from Pb-BSCCO (OD 72K) measured with 
21.2 eV photons with a low degree of Unear polarization, (a) 
raw MDM; (b) and (c) self-normalized MDM's using either 
the integrated intensity (b) or high binding energy intensity 
(c) as the denominator of the normalization function. 

FIG. 7. Contour plots of typical Pb-BSCCO FX EDMs for 
two temperatures: 30 K (top) and 300 K (bottom). The 
left panels show the experimental EDMs, from which a back- 
ground has been subtracted and the right panels show the 
results of the simulation. 

FIG. 8. (a) Simulation of the FX Ef-MDC for four differ- 
ent temperatures with = 20 meV and Rk = 30 meV (10~^ 
A-i). (b) Shift of the MDC maximum from the true kp as 
a function of temperature for different values of the energy 
resolution. 

FIG. 9. Simulations of (a) the gradient Iint{k) for the FX 
direction in BSCCO for four different temperatures. Ru, — 20 
meV and Rk = 30 meV (10~^ A"^) (see Eq. 6) and (b) the 
shift of the maximum in gradient Iint{k) from the true as 
a funtion of temperature for different momentum resolutions. 
For details see text. 
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